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The  dynamic  response  of  a  free-standing  plate  subjected  to  a  blast  wave  is  studied  numerically  to  investigate  the  effects  of 
fluid-structure  interaction  (FSI)  in  blast  wave  mitigation.  Previous  work  on  the  FSI  between  a  blast  wave  and  a  free-standing 
plate  (Kambouchev,  N.,  et  al.,  2006,  “Nonlinear  Compressibility  Effects  in  Fluid-Structure  Interaction  and  Their  Implications 
on  the  Air-Blast  Loading  of  Structures,”  J.  Appl.  Phys.,  100(6),  p.  063519)  has  assumed  a  constant  atmospheric  pressure  at  the 
back  of  the  plate  and  neglected  the  resistance  caused  by  the  shock  wave  formation  due  to  the  receding  motion  of  the  plate. 

This  paper  develops  an  FSI  model  that  includes  the  resistance  caused  by  the  shock  wave  formation  at  the  back  of  the  plate. 

The  numerical  results  show  that  the  resistance  to  the  plate  motion  is  especially  pronounced  for  a  light  plate,  and  as  a  result,  the 
previous  work  overpredicts  the  mitigation  effects  of  FSI.  Therefore,  the  effects  of  the  interaction  between  the  plate  and  the 
shock  wave  formation  at  the  back  of  the  plate  should  be  considered  in  blast  wave  mitigation. 
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Fluid  Structure  Interactions  for 
Blast  Wave  Mitigation 

The  dynamic  response  of  a  free-standing  plate  subjected  to  a  blast  wave  is  studied 
numerically  to  investigate  the  effects  of  fluid-structure  interaction  ( FSI )  in  blast  wave 
mitigation.  Previous  work  on  the  FSI  between  a  blast  wave  and  a  free-standing  plate 
(Kambouchev,  N.,  et  ah,  2006,  “Nonlinear  Compressibility  Effects  in  Fluid-Structure 
Interaction  and  Their  Implications  on  the  Air-Blast  Loading  of  Structures,  ”  J.  Appl. 
Phys.,  1 00(  6 ),  p.  063519)  has  assumed  a  constant  atmospheric  pressure  at  the  back  of  the 
plate  and  neglected  the  resistance  caused  by  the  shock  wave  formation  due  to  the  reced¬ 
ing  motion  of  the  plate.  This  paper  develops  an  FSI  model  that  includes  the  resistance 
caused  by  the  shock  wave  formation  at  the  back  of  the  plate.  The  numerical  results  show 
that  the  resistance  to  the  plate  motion  is  especially  pronounced  for  a  light  plate,  and  as 
a  result,  the  previous  work  overpredicts  the  mitigation  effects  of  FSI.  Therefore,  the 
effects  of  the  interaction  between  the  plate  and  the  shock  wave  formation  at  the  back  of 
the  plate  should  be  considered  in  blast  wave  mitigation.  [DOI:  10.1115/1.4002758] 
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1  Introduction 

Explosions  occur  in  military  conflicts,  as  well  as  in  various 
industrial  applications,  such  as  the  petrochemical,  chemical,  or 
nuclear  industries.  Underwater  shock  or  air  blast  waves  induced 
by  explosions  in  fluid  media  can  result  in  injuries  to  military  and 
civilian  personnel,  as  well  as  property  damage.  Hence,  our  under¬ 
standing  and  ability  to  correctly  model  complex  fluid  structure 
interaction  phenomena  has  important  consequences  for  mitigation 
of  infrastructural  damage  and  human  injury.  Understanding  the 
FSI  between  a  blast  wave  and  a  free-standing  plate  is  a  very 
important  step  in  the  development  of  devices  for  blast  wave  miti¬ 
gation.  A  blast  wave  impinging  on  a  free-standing  plate  will  cause 
the  plate  to  recede.  The  receding  motion  of  the  plate  relieves  the 
pressure  experienced  by  the  plate  and  results  in  a  decrease  in  the 
impulse  transmitted  to  the  plate.  Meanwhile,  the  plate  motion  in¬ 
duces  a  shock  wave  behind  the  plate,  which  resists  the  plate  mo¬ 
tion.  The  resistance  to  the  plate  motion  offsets  a  fraction  of  the 
decrease  in  pressure  experienced  by  the  plate,  which  is  due  to  the 
receding  motion  of  the  plate.  As  a  result,  the  impulse  transferred 
to  the  plate  increases  due  to  the  induced  shock  wave.  Both  the 
receding  motion  of  the  plate  and  the  resistance  to  the  plate  should 
be  considered  in  the  investigation  of  FSI  between  a  blast  wave 
and  a  free-standing  plate. 

In  recognition  of  the  advantages  of  FSI  in  blast  wave  mitiga¬ 
tion,  both  numerical  and  experimental  efforts  have  been  devoted 
to  this  research  area  during  the  past  decade.  Kambouchev  et  al. 
[1-3]  studied  the  FSI  between  a  blast  wave  and  a  free-standing 
plate  in  highly  compressible  medium,  having  assumed  that  a  con¬ 
stant  atmospheric  pressure  is  applied  at  the  back  of  the  plate.  The 
results  show  that  the  impulse  transmitted  to  structures  will  de¬ 
crease  with  the  reduction  of  the  plate  mass  for  fixed  blast  inten¬ 
sity.  Kambouchev  [4]  also  studied  the  interaction  of  acoustic 
waves  and  structures  in  incompressible  media.  The  acoustic  waves 
are  governed  by  linear  wave  equations.  The  effects  of  incompress¬ 
ible  media  on  both  sides  of  structures  were  considered,  and  an 
analytical  solution  was  derived.  Chun  et  al.  [5]  investigated  the 
FSI  of  flexible  shelters  under  blast  loading  using  ANSYS  and  FLU¬ 
ENT.  The  results  were  compared  with  experimental  data.  It  was 
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found  that  the  blast  loading  and  the  structural  response  cannot  be 
analyzed  separately  due  to  the  FSI  between  the  flexible  structure 
and  the  blast  loading.  Therefore,  FSI  should  be  considered  when 
the  structural  deflection  rate  can  influence  the  solution  of  the  flow 
field  surrounding  the  structure.  Main  and  Gazonas  [6]  studied  the 
influence  of  mass  distribution  on  the  uniaxial  crushing  of  a  cellu¬ 
lar  material  sandwiched  between  solid  front  and  back  faces,  with 
air  blast  loading  applied  to  the  front  face  and  the  back  face  unre¬ 
strained.  An  analytical  model  was  developed  to  investigate  the 
crushing  response  of  the  system.  FSI  effects  were  treated  using  the 
results  of  Kambouchev  et  al.  [1,2]  for  air  blast  loading  on  solid 
plates.  It  was  found  that  the  FSI  approximation  of  Kambouchev  et 
al.  [1,2]  is  inappropriate  when  the  mass  fraction  in  the  core  is 
large  relative  to  the  mass  in  the  front  face,  particularly  when  the 
decay  period  of  the  pressure  pulse  is  comparable  to  the  propaga¬ 
tion  time  of  the  densification  front.  Espinosa  et  al.  [7]  proposed  a 
novel  experimental  methodology  that  incorporates  FSI  to  assess 
the  dynamic  deformation  of  structures  subjected  to  underwater 
impulsive  loading. 

Related  papers  studying  the  blast  wave  propagation  and  reflec¬ 
tion  have  also  appeared  in  the  literature.  Tai  et  al.  [8]  discussed 
blast  wave  interaction  and  reflection  around  closed-ended  and 
open-ended  bomb  shelters.  The  reflected  shock  wave  patterns 
transit  from  regular  reflection  to  Mach  reflection  in  both  bomb 
shelters  under  an  unsteady  situation.  Shi  et  al.  [9]  simulated  the 
blast  wave  interaction  with  a  standalone  structural  column.  Para¬ 
metric  studies  include  the  scaled  distance  of  the  blast,  column 
stiffness,  and  column  dimension  and  geometry.  The  formula  to 
predict  the  reflected  pressure  and  impulse  on  the  front  and  on  the 
rear  surface  of  the  columns  with  different  dimensions  and  geom¬ 
etry  were  derived.  Blast  wave  propagation  and  reflection  has  been 
investigated  experimentally.  Dewey  and  McMillin  [10]  used  high 
speed  photogrammetry  to  investigate  the  blast  wave  interaction 
with  ideal  and  real  surfaces.  It  was  observed  that  a  smooth  surface 
induces  a  stronger  Mach  stem.  Takayama  and  Sekiquchi  [11]  stud¬ 
ied  the  interaction  of  a  spherical  shock  wave  with  a  planar  or 
conical  wall.  The  shock  wave  was  induced  by  expanding  a  planar 
shock  wave  into  free  space  in  a  conventional  shock  tube.  When  a 
spherical  shock  wave  encounters  a  planar  or  conical  wall,  a  tran¬ 
sition  from  regular  to  Mach  reflection  takes  place  with  the  inci¬ 
dence  angle  larger  than  the  critical  transition  angle.  Igra  et  al.  [12] 
investigated  blast  wave  reflection  from  wedges  experimentally 
and  numerically.  A  shock  tube  equipped  with  a  very  short  driver  is 
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Fig.  1  Schematic  of  a  free-standing  plate  subjected  to  a  blast 
wave 


employed  for  the  experiments.  The  reflected  wave  pattern  is  simi¬ 
lar  for  the  interaction  of  a  blast  wave  or  a  shock  wave  with  a 
wedge  when  both  incident  waves  have  the  same  initial  pressure 
jump  across  their  fronts,  but  the  resulting  pressure  field  is  differ¬ 
ent. 

Previous  work  [1-3]  on  the  FSI  between  a  blast  wave  and  a 
free-standing  plate  in  highly  compressible  medium  has  neglected 
the  shock  wave  induced  at  the  back  of  the  plate  caused  by  plate 
motion.  This  assumption  leads  to  overprediction  of  the  effective¬ 
ness  of  FSI  in  blast  wave  mitigation.  In  another  work  [4],  the 
effects  of  incompressible  media  on  both  sides  of  structures  were 
considered  for  FSI  between  acoustic  waves  and  structures.  Since 
the  blast  wave  is  a  shock  wave  governed  by  nonlinear  wave  equa¬ 
tions  as  opposed  to  an  acoustic  wave  governed  by  linear  wave 
equations,  their  simplification  results  in  conclusions  that  cannot  be 
readily  applied  to  blast  wave  mitigation.  To  accurately  model  the 
FSI  between  a  blast  wave  and  a  free-standing  plate  in  highly  com¬ 
pressible  medium,  it  is  necessary  to  develop  a  model  that  includes 
the  shock  wave  induced  by  the  plate  motion.  In  this  paper,  a 
numerical  model  is  developed  to  simulate  the  blast  wave  reflec¬ 
tion  from  a  free-standing  plate,  the  plate  motion,  and  the  shock 
wave  induced  at  the  back  of  the  plate.  The  model  is  validated  by 
comparing  results  with  numerical  and  analytical  results  available 
in  the  literature  for  special  cases.  The  numerical  results  are  pre¬ 
sented  and  analyzed  in  subsequent  sections. 


2  Mathematical  Formulation  and  Numerical  Ap¬ 
proach 


2.1  Governing  Equations.  The  schematic  of  the  FSI  between 
a  blast  wave  and  a  free-standing  plate  is  shown  in  Fig.  1 .  The  plate 
separates  the  computational  domain  into  two  parts:  flow  field  in 
front  of  the  plate  and  flow  field  at  the  back  of  the  plate.  The  flow 
fields  on  both  sides  can  be  modeled  as  one-dimensional  inviscid 
compressible  flows,  which  are  described  by  the  Euler  equations. 
They  can  be  written  in  vector  form  as 


dU  dF 
—  +  —  =  0 
at  ax 


(i) 


where  U  is  the  solution  vector,  F  is  the  flux  vector,  t  is  time,  and 
x  is  the  space  coordinate  in  the  flow  direction.  U  and  F  are  given 
by 


p 

pu 

u= ■ 

pu 

'  F  = 

pu 1  +  p 

pe, 

( pe,+p)u 

where  p  is  the  density,  p  is  the  pressure,  u  is  the  velocity,  and  e,  is 
the  total  energy  per  unit  mass  of  the  compressible  flow.  The  total 
energy  per  unit  mass  is  the  sum  of  its  internal  energy  per  unit 
mass,  e,  and  its  kinetic  energy  per  unit  mass,  irl  2.  To  close  the 
governing  equations,  it  is  assumed  that  the  compressible  sub¬ 
stance  obeys  calorically  perfect  ideal  gas  law  given  by  p=pRT 
and  e  =  (RT)l (y-1),  where  T  is  the  temperature,  R  is  the  gas  con¬ 


stant,  and  y  is  the  specific  heat  ratio. 

The  free-standing  plate  is  treated  as  a  rigid  body,  and  the  effects 
of  deformation  and  stress-wave  propagation  within  the  plate  are 
neglected.  The  free-standing  plate  obeys  Newton’s  second  law  of 
motion 


dup  A pp 
dt  pphp 


(2) 


where  up  is  the  velocity  of  the  plate,  pp  is  the  density  of  the  plate, 
hp  is  the  thickness  of  the  plate,  and  A pp  is  the  difference  between 
the  reflected  pressure  in  front  of  the  plate  and  the  induced  shock 
wave  pressure  at  the  back  of  the  plate. 


2.2  Initial  and  Boundary  Conditions.  The  initial  conditions 
for  the  flow  fields  (p0,T0,u0)  are  prescribed  as  the  ambient  con¬ 
ditions.  The  free-standing  plate  is  initially  at  rest  {up= 0  m/s). 

Two  types  of  blast  waves  are  studied  in  this  paper.  The  uniform 
incident  blast  wave  can  be  written  as 


Pi(t)  =  Pu  0) 

where  pu  is  a  constant  overpressure  of  the  incident  blast  wave. 
The  second  type  simulates  a  typical  exponential  blast  wave,  which 
consists  of  an  abrupt  pressure  increase  followed  by  slow  decay  of 
pressure.  There  is  usually  a  minor  negative  phase  at  the  tail  end  of 
the  blast  wave.  Neglecting  the  negative  phase,  the  exponential 
incident  blast  wave  can  be  approximated  by  an  exponential  profile 

Pi(t)  =  Pise~,"i  (4) 

where  pis  is  the  peak  overpressure  of  the  incident  blast  wave  that 
is  at  a  stand-off  distance  rf,  from  the  free-standing  plate  and  r,-  is 
the  decay  time  constant.  This  exponential  profile  does  not  corre¬ 
spond  to  real  blast  wave  propagation  but  to  an  approximation  that 
is  used  for  initialization.  The  density  and  velocity  of  both  types  of 
incident  blast  waves  are  related  to  the  pressure  through  the  clas¬ 
sical  Rankine-Hugoniot  relations 


Pi(t )  =  Po 


(y+  l)pj(f)  +  2yp0 

(y-  i)pM  +  2ypo 


«,-(0 : 


a0Pi(t) 


yPo  v(y+  l)pj(t)  +  2yp0 


(5) 


where  a0  is  the  speed  of  sound  at  ambient  conditions. 

The  above  equations  serve  as  the  left-hand  side  boundary  con¬ 
ditions  for  the  flow  field  in  front  of  the  plate.  On  the  right-hand 
side  of  this  flow  field,  zero  gradient  boundary  conditions  are  ap¬ 
plied  to  pressure  and  density.  The  no-penetration  condition  is  ap¬ 
plied  to  the  velocity.  These  boundary  conditions  can  be  written  as 


u  =  up  (6) 

The  computational  domain  for  the  flow  field  in  front  of  the  plate 
spans  from  the  initial  blast  wave  to  the  free-standing  plate. 

The  left-hand  side  boundary  conditions  for  the  flow  field  at  the 
back  of  the  plate  are  the  same  as  the  right-hand  side  boundary 
conditions  for  the  flow  field  in  front  of  the  plate.  The  right-hand 
side  boundary  conditions  for  the  flow  field  at  the  back  of  the  plate 
are  given  by  the  ambient  conditions  and  can  be  written  as 


P=Po 


P  =  Po 
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U  =  Uq 


(7) 

The  initial  conditions  for  the  flow  fields  on  both  sides  of  the  plate 
can  be  described  by  Eq.  (7).  The  computational  domain  for  the 
flow  field  at  the  back  of  the  plate  should  be  large  enough  so  that 
the  right  boundary  conditions  will  not  be  affected  by  the  induced 
shock  wave. 

2.3  Numerical  Approach.  The  Euler  equations  for  the  flow 
fields  both  in  front  and  at  the  back  of  the  plate  are  solved  using  the 
van  Leer  flux  vector  splitting  scheme  [13,14]  coupled  with  the 
monotone  upstream-centered  scheme  for  conservation  laws 
(MUSCL)  and  Runge-Kutta  scheme.  Because  there  are  sharp  gra¬ 
dients  in  the  flow  fields  for  blast  wave  reflection  and  shock  wave 
generation,  numerical  dissipation  is  needed  in  the  numerical 
scheme  to  attenuate  numerical  errors  of  small  wavelengths.  A 
standard  upwind  scheme  can  provide  sufficient  numerical  dissipa¬ 
tion  for  most  flow  problems  that  can  involve  shock  waves  with 
large  gradients.  However,  the  directions  of  characteristic  velocity 
should  be  established  before  the  upwind  scheme  can  be  imple¬ 
mented.  The  van  Leer  flux  vector  splitting  scheme  is  applied  to 
split  the  flux  vector  F  into  forward  and  backward  components  F+ 
and  F~  in  terms  of  the  local  Mach  number.  F+  and  F~  have  char¬ 
acteristic  velocities  in  forward  and  backward  directions,  respec¬ 
tively.  MUSCL  is  then  used  to  discretize  the  spatial  derivative  of 
the  flux  vector.  MUSCL  is  a  first-order  upwind  scheme  in  the 
vicinity  of  the  shock  and  a  second-order  upwind  scheme  else¬ 
where  in  the  flow  field.  It  provides  good  stability  as  well  as  accu¬ 
racy.  A  second-order  Runge-Kutta  explicit  scheme  is  used  to  dis¬ 
cretize  the  temporal  derivatives  of  the  Euler  equations. 

A  second-order  Runge-Kutta  explicit  scheme  is  used  to  dis¬ 
cretize  Newton’s  second  law  of  motion  for  the  free-standing  plate. 
The  detailed  numerical  scheme  and  solution  algorithm  have  been 
described  in  an  earlier  work  [15]. 


3  Numerical  Validation 

To  verify  the  accuracy  of  the  numerical  model  developed  in  this 
study,  the  flow  fields  in  front  and  at  the  back  of  the  plate  were 
simulated  separately  and  compared  with  results  available  in  the 
literature. 

The  following  parameters  were  used  in  the  numerical  simula¬ 
tion.  The  specific  heat  ratio  y  and  the  gas  constant  R  for  air  are 
1.4  J/kg  K  and  287  J/kg  K,  respectively.  The  density  of  the  steel 
plate  pp  is  7800  kg/m3.  The  ambient  conditions  are 

p0  =  101.3  kPa 
T0  =  298  K 
u0  =  0  m/s 

3.1  Model  Validation:  Part  I.  The  simulation  of  the  com¬ 
pressible  flow  field  in  front  of  the  free-standing  plate,  which  is 
subjected  to  an  exponential  blast  wave,  is  presented  first.  The 
resistance  to  the  plate  motion  is  neglected  by  assuming  a  constant 
atmospheric  pressure  at  the  back  of  the  plate.  The  FSI  problem 
under  this  assumption  was  investigated  by  Kambouchev  et  al.  [1], 
Based  on  the  acoustic  results  for  low-intensity  waves  and  the 
asymptotic  limits  for  extremely  light  and  extremely  heavy  plates 
placed  within  an  air  environment,  they  developed  a  semi- 
analytical  formula  to  predict  the  ratio  of  impulse  transmitted  to  the 
plate  (Ip)  to  the  incident  impulse  (/,■)  for  arbitrary  blast  intensity 
and  plate  mass.  The  impulse  transmitted  to  the  plate,  defined  as 
Ip=S(Pr-Pb)dt,  was  used  by  Kambouchev  et  al.  [1]  and  shows  the 
net  impulse  transmitted  to  the  plate.  Here,  pr  is  the  overpressure 
of  the  reflected  blast  wave  in  front  of  the  plate  and  pb  is  the 
overpressure  at  the  back  of  the  plate.  The  incident  impulse  is 
defined  as  Ii=Jpldt,  where  p,  is  the  overpressure  of  the  incident 
blast  wave. 


Fig.  2  The  ratio  of  impulse  versus  compressible  FSI  param¬ 
eter  for  different  blast  intensities 
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Here,  pis  is  the  peak  density  behind  the  incident  shock  front,  Uis  is 
the  peak  propagation  speed  of  the  incident  blast  wave,  and  f3s  is  a 
nondimensional  parameter  describing  FSI  and  is  inversely  propor¬ 
tional  to  the  plate  thickness  hp. 

The  numerical  results  in  this  paper  are  compared  with  the  re¬ 
sults  of  Kambouchev  et  al.  [1]  for  two  different  blast  intensities  of 
3.29  and  10.85  and  at  six  different  plate  thicknesses,  as  shown  in 
Fig.  2.  It  can  be  seen  that  the  numerical  results  of  the  present 
study  agree  well  with  the  previous  results  of  Kambouchev  et  al. 
[1].  The  ratio  of  impulse  I p / 1,  decreases  with  an  increase  in  /3S 
(decrease  in  plate  thickness).  For  a  plate  with  large  thickness 
(small  f)s ),  it  behaves  like  a  fixed  wall.  The  corresponding  ratio  of 
impulse  lpH,  approaches  that  of  a  fixed  wall.  The  reduction  in  the 
ratio  of  impulse  I p /  /,  is  minimal.  For  relatively  thick  plates,  there 
is  no  significant  blast  wave  mitigation.  For  thin  plates  (large  /3S), 
the  ratio  of  impulse  IpIIj  becomes  very  small  for  both  blast  inten¬ 
sities.  This  means  that  the  blast  wave  mitigation  using  a  thin  plate 
is  very  effective. 


3.2  Model  Validation:  Part  II.  The  flow  field  behind  the 
plate  is  simulated  by  including  the  resistance  to  the  plate  motion. 
A  shock  wave  is  induced  at  the  back  of  the  plate  due  to  the  plate 
motion.  When  the  free-standing  plate  in  the  left  boundary  moves 
with  a  constant  velocity  up,  the  overpressure  of  the  shock  wave 
induced  by  the  plate  motion  is  available  analytically  [16,17], 
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Fig.  3  Numerical  and  analytical  pressure  profiles  of  the  induced  shock  wave 
behind  the  plate  for  (a)  low  and  ( b )  high  plate  velocities  at  different  times 
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where  pg  is  the  overpressure  of  the  shock  wave  induced  by  the 
plate  motion.  With  a  known  value  of  up,  the  overpressure  pg  can 
be  determined  using  Eq.  (13). 

A  plate  moving  at  a  constant  velocity  induces  a  shock  wave, 
which  propagates  at  a  higher  speed  in  the  same  direction  as  the 
plate  motion.  The  pressure  profiles  of  the  induced  shock  wave 
obtained  from  the  FSI  model  are  compared  with  the  analytical 
solution  at  different  times  of  0.1  ms  and  0.5  ms  and  for  plate 
velocities  of  100  m/s  and  10,000  m/s.  These  two  plate  velocities 
are  used  to  demonstrate  the  validity  of  the  numerical  model.  It  can 
be  observed  that  the  numerical  results  are  in  good  agreement  with 
the  analytical  solution.  The  starting  point  of  each  curve  at  differ¬ 
ent  times  of  0.1  ms  and  0.5  ms  indicates  the  position  of  the  free¬ 
standing  plate.  As  shown  in  Fig.  3,  the  shock  wave  strength  in¬ 
duced  by  the  plate  velocity  of  10,000  m/s  is  many  orders  of 
magnitude  larger  than  that  induced  by  the  plate  velocity  of  100 
m/s. 


3.3  Model  Validation:  Part  III.  The  flow  fields  on  both  sides 
of  the  plate  are  coupled  through  the  receding  motion  of  the  plate, 
which  induces  a  shock  wave  at  the  back  of  the  plate.  In  Ref.  [4], 
Kambouchev  derived  an  analytical  solution  for  acoustic  waves 
interacting  with  free-standing  plates  with  the  presence  of  a  sup¬ 
porting  acoustic  medium  at  the  back  of  the  plate.  The  solution  is 
applicable  for  acoustic  waves  that  are  governed  by  linear  wave 
equations.  The  net  impulse  transmitted  to  the  plate  is  given  by  [4] 

/  ft 

f  =  2/0‘-ft  (14) 

where  /30  is  the  fluid-structure  interaction  parameter  given  by 
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Here,  a  is  the  speed  of  sound  in  the  medium  in  front  of  the  plate 
whose  density  is  p,  ar  is  the  speed  of  sound  in  the  medium  at  the 
back  of  the  plate  whose  density  is  ~pr,  f,  is  the  decay  time  of  the 
exponential  incident  blast  wave,  and  mp  is  the  mass  of  the  plate 
per  unit  area. 

At  low  blast  intensity,  the  shock  wave  behaves  as  an  acoustic 
wave.  Numerical  results  obtained  from  the  present  model  under 
such  condition  can  be  compared  with  the  analytical  solution  of 
Kambouchev  [4],  Numerical  simulations  have  been  carried  out  for 
six  different  plate  thicknesses  at  an  incident  acoustic  wave,  which 
has  a  peak  intensity  of  0.2  and  a  decay  time  of  1  ms.  The  com¬ 
parison  between  the  numerical  results  and  the  analytical  solution 


of  Kambouchev  [4]  is  shown  in  Fig.  4.  It  can  be  seen  that  the 
numerical  results  agree  well  with  Kambouchev’s  analytical  solu¬ 
tion. 

4  Results  and  Discussion 

Numerical  simulations  of  the  FSI  between  a  blast  wave  and  a 
free-standing  plate  have  been  carried  out.  The  shock  wave  in¬ 
duced  at  the  back  of  the  plate  is  taken  into  account.  In  this  model, 
the  flow  fields  on  both  sides  of  the  plate  are  coupled  through  the 
receding  motion  of  the  plate. 

4.1  Uniform  Blast  Wave.  As  discussed  before,  the  FSI  be¬ 
tween  the  receding  plate  and  the  induced  shock  wave  at  the  back 
of  the  plate  affects  the  blast  wave  reflection.  The  reflection  coef¬ 
ficient  C  is  typically  defined  as  the  ratio  of  the  reflected  blast 
wave  overpressure  to  the  incident  blast  wave  overpressure.  The 
reflection  coefficient  is  analyzed  for  a  free-standing  plate  sub¬ 
jected  to  a  uniform  blast  wave.  There  are  two  limiting  cases  of  this 
physical  problem.  One  corresponds  to  a  plate  of  infinite  mass,  and 
the  other  one  corresponds  to  a  plate  of  infinitesimal  mass.  For  a 
free-standing  plate  of  infinite  mass,  the  reflection  coefficient  cor¬ 
responds  to  that  for  a  fixed  wall,  which  has  been  derived  previ¬ 
ously  and  can  be  written  as  [16] 

c_  Pr  _  (3y~  l)pu/Po  +  4y  ^ 

Pu  (y-i)pJpo+2y 

For  air,  the  range  of  the  reflection  coefficient  for  a  plate  of  infinite 
mass  is 


2<C<8 

The  minimum  value  of  2  corresponds  to  an  acoustic  reflection, 


Fig.  4  The  ratio  of  impulse  versus  FSI  parameter  for  acoustic 
wave 
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and  the  maximum  value  of  8  corresponds  to  a  strong  blast  wave 
reflection. 

For  a  free-standing  plate  of  infinitesimal  mass,  the  inertial  ef¬ 
fect  of  the  plate  is  negligible.  As  the  blast  wave  impacts  the  plate, 
the  plate  accelerates  instantaneously  to  the  equilibrium  state.  This 
implies  that  the  force  exerted  by  the  reflected  blast  wave  in  front 
of  the  plate  should  be  the  same  as  that  exerted  by  the  induced 
shock  wave  at  the  back  of  the  plate. 

As  discussed  in  Sec.  3.2,  the  shock  wave  strength  pglpa  in¬ 
duced  by  the  plate  motion  is  related  to  the  plate  velocity  up 
through  Eq.  (13).  The  reflected  blast  wave  strength  is  related  to 
the  plate  velocity  up  through  the  Rankine-Hugoniot  relations.  It 
should  be  noted  that  the  Rankine-Hugoniot  relations  for  a  uni¬ 
form  reflected  blast  wave  are  different  from  the  Rankine- 
Hugoniot  relations  for  a  uniform  incident  blast  wave.  They  can  be 
expressed  as 


where  uu 
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is  the  fluid  velocity  behind  the  incident  shock  front 
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and  Wr  is  the  reflected  blast  wave  speed  given  by 

Wr  =  arMr-uu  (19) 

The  Mach  number  of  the  reflected  blast  wave,  Mn  is  given  by 


M,.= 


(p_lZPu 

\Pu+P0 


+  1 


(20) 


and  the  speed  of  sound  ar  is  related  to  Tln  the  temperature  behind 
the  incident  shock  front,  through  the  equation 


T„  =  T, 
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For  a  free-standing  plate  of  infinitesimal  mass,  pg  is  equal  to  pr. 
Combining  Eqs.  (17)— (21)  and  Eq.  (13)  yields  a  reflection  coeffi¬ 
cient  C  of  1. 

A  free-standing  plate  initially  at  rest  behaves  like  a  fixed  wall. 
The  instant  that  the  blast  wave  impacts  the  plate,  the  reflection 
coefficient  C  is  expected  to  approach  the  value  corresponding  to  a 
fixed  wall.  The  reflection  coefficient  C  then  gradually  decreases  to 
its  equilibrium  value.  Numerical  simulations  have  been  carried 
out  for  free-standing  plates  of  5  mm  and  50  mm  thickness  at 
uniform  blast  intensities  of  5  and  50.  Figure  5  shows  how  the 
reflection  coefficient  C  varies  with  time.  For  a  fixed  wall,  the 
reflection  coefficients  C  for  uniform  blast  wave  with  blast  inten¬ 
sities  of  5  and  50  are  4.5  and  7.2,  respectively.  As  shown  in  Fig.  5, 
at  the  instant  that  the  blast  wave  impacts  the  plate,  the  reflection 
coefficients  of  free-standing  plates  approach  4.5  and  7.2  for  uni¬ 
form  blast  intensities  of  5  and  50,  respectively.  The  reflection 
coefficient  gradually  decreases  with  time  to  its  equilibrium  value. 
The  equilibrium  value  for  the  thinner  plate  is  lower  than  that  for 
the  thicker  plate  due  to  the  FSI  between  the  blast  wave  and  the 
plate. 

To  study  how  the  equilibrium  value  of  reflection  coefficient  C 
varies  with  uniform  blast  intensities,  numerical  simulations  of  the 
FSI  have  been  carried  out  for  four  different  plate  thicknesses: 
infinite  thickness,  50  mm,  15  mm,  and  5  mm.  Figure  6  shows  the 


Fig.  5  The  reflection  coefficient  prlpu  as  a  function  of  time  t 


equilibrium  value  of  reflection  coefficient  C  for  all  simulated  plate 
thicknesses  at  different  uniform  incident  blast  intensities.  In  Fig. 
6,  the  numerical  results  for  a  plate  of  infinite  mass  are  compared 
with  the  analytical  solution  that  corresponds  to  this  limiting  case. 
It  can  be  seen  that  the  numerical  results  for  an  infinite  plate  thick¬ 
ness  agree  well  with  the  analytical  results  for  a  fixed  wall.  The 
reflection  coefficient  increases  from  2  to  8  with  the  increase  of 
incident  blast  intensity.  For  a  small  incident  blast  intensity,  all 
free-standing  plates  are  relatively  heavy  to  the  incident  blast 
wave.  Therefore,  the  reflection  coefficient  for  all  free-standing 
plates  approaches  asymptotically  that  of  a  fixed  wall.  Figure  6 
also  shows  that  as  the  incident  blast  intensity  increases,  the  reflec¬ 
tion  coefficient  deviates  from  that  of  the  fixed  wall  gradually  due 
to  the  FSI.  The  inception  of  the  deviation  for  a  thinner  plate  starts 
earlier  than  that  for  a  thicker  plate.  There  is  a  critical  transition 
point  for  a  particular  plate  thickness.  When  the  incident  blast  in¬ 
tensity  is  greater  than  the  critical  transition  value,  the  reflection 
coefficient  starts  to  decrease.  It  indicates  that  the  FSI  plays  a  more 
important  role  for  large  blast  intensities.  As  the  blast  intensity 
increases  further,  the  reflection  coefficient  of  the  free-standing 
plate  of  finite  mass  gradually  approaches  to  that  of  a  plate  of 
infinitesimal  mass.  As  shown  in  Fig.  6(a),  the  reflection  coefficient 
for  the  5  mm  thick  plate  approaches  1  asymptotically  when  the 
blast  intensity  exceeds  30.  The  reflection  coefficients  for  plates  of 
15  mm  and  50  mm  approach  1  as  the  blast  intensity  continues  to 
increase  (see  Fig.  6(b)).  It  should  be  noted  that  the  reflection  co¬ 
efficient  C  can  never  be  less  than  1. 

4.2  Exponential  Blast  Wave.  Numerical  simulations  of  an 
exponential  blast  wave  interacting  with  a  free-standing  plate  have 
been  conducted  for  different  blast  intensities  and  plate  thick¬ 
nesses.  Because  both  the  incident  pressure  and  reflected  pressure 
vary  with  time  for  the  exponential  blast  wave,  the  reflection  coef¬ 
ficient  as  defined  in  Sec.  4.1  is  not  directly  applicable  to  the  ex¬ 
ponential  blast  wave.  Therefore,  the  reflected  overpressure  pr  and 
the  incident  overpressure  p;  are  integrated  over  the  duration  of  the 
blast  wave  to  obtain  the  impulse  transmitted  to  the  plate  Ip 
=  fprdt  and  the  incident  impulse  Ij=fpjdt.  Physically,  the  ratio  of 
impulse  Ipll,  is  the  average  of  the  reflection  coefficient  over  time. 

As  can  be  seen  in  the  previous  paragraph,  the  transmitted  im¬ 
pulse  used  in  this  paper  is  different  from  the  definition  used  by 
Kambouchev  et  al.  [1],  which  considered  the  net  impulse  trans¬ 
mitted  to  the  plate,  namely,  Ip=f(pr-pb)dt.  The  definition  in  the 
present  study  considers  the  impulse  transmitted  to  the  plate  by  the 
blast  wave  and  properly  represents  the  damage  potential  of  the 
blast  wave.  This  definition  is  more  meaningful  in  analysis  of  blast 
wave  impact.  It  can  be  illustrated  by  a  simple  thought  experiment. 
For  example,  if  a  soldier  is  impacted  by  two  comparable  blast 
waves  on  the  front  and  at  the  back  simultaneously,  the  soldier  will 
suffer  from  grave  injury  or  death.  The  net  impulse  transmitted  to 
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Fig.  6  The  reflection  coefficient  for  different  plate  thicknesses  at  different  uni¬ 
form  incident  blast  intensities,  (a)  Uniform  incident  blast  intensities  of  0-50;  (b) 
uniform  incident  blast  intensities  of  0-500. 


the  soldier,  using  the  definition  of  Kambouchev  et  al.  [1],  is  neg¬ 
ligible.  The  net  transmitted  impulse  does  not  provide  proper  mea¬ 
surement  of  the  damage  caused  by  the  blast  wave  because  the 
damage  potential  of  two  blast  waves  cannot  be  cancelled.  On  the 
other  hand,  using  the  definition  in  this  study,  both  impulses  trans¬ 
mitted  to  the  soldier  by  the  blast  waves  are  significant  and  provide 
a  good  measurement  of  the  damage  caused  by  the  blast  waves. 

A  blast  wave  with  a  peak  overpressure  of  2  MPa,  decay  time 
constant  of  0.5  ms,  and  stand-off  distance  of  1.6  m  is  adopted  for 
the  numerical  simulation.  This  results  in  a  peak  blast  intensity  of 
10.85  right  in  front  of  the  plate.  Approximately  the  same  plate 
thicknesses  (0.000378  mm,  0.00378  mm,  0.0378  mm,  0.378  mm, 
3.78  mm,  and  37.8  mm)  as  in  Ref.  [1]  have  been  simulated.  The 
ratio  of  impulse  transmitted  to  the  plate  ( Ip )  to  the  incident  im¬ 
pulse  (/,■)  for  low  blast  intensity  of  10.85  with  resistance  at  the 
back  of  the  plate  is  shown  in  Fig.  7.  For  relatively  heavy  plates, 
the  impact  of  incident  blast  wave  results  in  a  relatively  small 
receding  velocity.  This  implies  that  the  blast  wave  mitigation  ef¬ 
fects  of  FSI  for  relatively  heavy  plates  are  insignificant.  As  the 
thickness  of  the  plate  decreases,  the  receding  velocity  of  the  plate 
increases.  The  ratio  of  impulse  IpIIj  decreases.  For  thicknesses 
less  than  about  0.378  mm,  the  ratio  of  impulse  Ipl Ij  reaches  a 
plateau.  This  is  due  to  the  increased  resistance  at  the  back  of  the 
plate.  The  impact  of  the  blast  wave  drives  the  relatively  light  plate 
to  very  high  receding  velocity,  which  induces  a  strong  shock  wave 
at  the  back  of  the  plate.  When  the  shock  wave  strength  at  the  back 
of  the  plate  becomes  comparable  to  the  incident  blast  intensity,  the 
effectiveness  of  FSI  in  reducing  the  blast  wave  impulse 


Fig.  7  Comparison  between  the  ratio  of  impulse  /p//,  with  and 
without  resistance  for  low  blast  intensity  at  different  plate 
thicknesses 


diminishes. 

The  ratio  of  impulse  I p /  /,  was  investigated  in  Ref.  [1],  where  Ip 
was  defined  as  Ip  =  J(pr-pb)dt.  Since  the  resistance  to  the  plate 
motion  was  neglected  by  assuming  a  constant  atmospheric  pres¬ 
sure  at  the  back  of  the  plate  in  the  previous  model  [1],  the  over¬ 
pressure  at  the  back  of  the  plate  pb  is  0.  Therefore,  Ip=f(pr 
-pb)dt  equals  Ip=Jprdt.  From  the  mathematical  point  of  view,  the 
ratio  of  impulse  I p  /  /,  defined  in  this  paper  is  the  same  as  that  in 
Ref.  [1],  To  investigate  whether  the  resistance  at  the  back  of  the 
plate  is  important,  the  difference  between  the  numerical  results 
obtained  from  the  two  different  models  is  studied.  The  FSI  be¬ 
tween  a  blast  wave  and  a  free-standing  plate  without  considering 
shock  wave  formation  at  the  back  of  the  plate  [1]  is  also  shown  in 
Fig.  7.  For  relatively  heavy  plates,  the  effects  of  resistance  on  FSI 
are  not  very  important.  The  difference  between  the  ratio  of  im¬ 
pulse  /p/7;  with  and  without  resistance  is  very  small.  As  the  thick¬ 
ness  of  the  plate  decreases,  the  effects  of  resistance  increase.  The 
difference  between  the  ratio  of  impulse  Jp!  1,  with  and  without 
resistance  becomes  significant.  While  the  previous  model  [1]  pre¬ 
dicts  that  it  is  possible  to  achieve  100%  reduction  with  a  light 
free-standing  plate  by  neglecting  the  resistance  at  the  back  of  the 
plate,  the  ratio  of  impulse  Ipl It  is  limited  by  the  resistance  at  the 
back  of  the  plate  to  around  2  for  a  blast  intensity  of  10.85.  This 
agrees  with  the  results  obtained  in  Sec.  4.1.  It  can  be  seen  from 
Fig.  6  that  the  reflection  coefficient  of  a  free-standing  plate  can 
never  be  less  than  1.  Since  the  ratio  of  impulse  transmitted  to  the 
plate  to  the  incident  impulse  /p/7,  is  the  average  reflection  coeffi¬ 
cient  C,  the  ratio  of  impulse  lpUj  is  expected  to  be  always  larger 
than  1.  Hence,  the  previous  model  [1]  significantly  overpredicts 
the  effectiveness  of  FSI  in  reducing  the  blast  wave  impact.  For 
this  low  blast  intensity,  the  range  of  small  plate  thicknesses  where 
the  two  physical  models  differ  in  blast  wave  mitigation  is  of  lim¬ 
ited  practical  value.  However,  the  difference  becomes  important 
for  more  realistic  conditions  presented  in  Figs.  8  and  9. 

Results  for  a  blast  intensity  of  100  are  presented  in  Fig.  8.  A 
blast  wave  with  a  peak  overpressure  of  16  MPa,  decay  time  con¬ 
stant  of  0.5  ms,  and  stand-off  distance  of  3.2  m  is  adopted.  These 
parameters  result  in  a  peak  blast  intensity  of  100  right  in  front  of 
the  plate.  Five  different  plate  thicknesses  ranging  from  0.00378 
mm  to  37.8  mm  have  been  simulated.  The  FSI  between  a  blast 
wave  and  a  free-standing  plate  for  a  blast  intensity  of  100  behaves 
very  much  similar  to  the  blast  intensity  of  10.85.  It  is  found  that 
the  ratio  of  impulse  /p//,  reaches  a  limit  of  3.6  for  thin  plates.  The 
effectiveness  of  FSI  in  blast  wave  mitigation  decreases  with  the 
increased  blast  intensity. 

Additional  simulations  are  conducted  with  a  blast  intensity  of 
1000.  This  is  achieved  by  placing  a  blast  wave  with  a  peak  over¬ 
pressure  of  150  MPa  and  decay  time  constant  of  0.5  ms  at  a 
stand-off  distance  of  9.6  m.  Four  different  plate  thicknesses  rang- 
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Fig.  8  Comparison  between  the  ratio  of  impulse  /p//,  with  and 
without  resistance  for  higher  blast  intensity  at  different  plate 
thicknesses 


Fig.  10  The  ratio  of  impulse  lpl /,  as  a  function  of  incident  im¬ 
pulse  /, 


ing  from  0.0378  mm  to  37.8  mm  have  been  simulated.  For  this 
blast  wave,  the  ratio  of  impulse  I pl I,  reaches  a  limit  of  5.4  for  thin 
plates.  The  difference  between  the  ratio  of  impulse  I p / /,  with  and 
without  resistance  becomes  100%  for  plate  thickness  of  3.78  mm. 
As  the  blast  intensity  reaches  1000,  the  relatively  thick  plate  of 
3.78  mm  is  relatively  light  for  the  blast  wave.  Therefore,  the  re¬ 
sistance  at  the  back  of  the  plate  significantly  affects  the  FSI.  Ne¬ 
glecting  the  resistance  causes  large  error  in  predicting  the  ratio  of 
impulse  lpHr 

The  ratio  of  impulse  Ip/ 1,  is  not  only  affected  by  the  incident 
peak  blast  intensity  but  also  by  the  incident  impulse.  For  a  par¬ 
ticular  incident  peak  blast  intensity,  different  decay  time  constants 
generate  different  incident  impulses.  Simulations  of  a  blast  wave 
interacting  with  a  free-standing  plate  of  3.78  mm  are  conducted 
for  different  incident  impulses.  An  incident  blast  wave  with  a  peak 
overpressure  of  8  MPa  and  decay  time  constants  ranging  from  0. 1 
ms  to  1.0  ms  is  placed  at  different  stand-off  distances  to  obtain  a 
constant  peak  blast  intensity  of  50  and  decay  time  constants  rang¬ 
ing  from  0.14  ms  to  1.9  ms  right  in  front  of  the  plate  with  different 
incident  impulses.  For  an  infinitesimal  impulse,  the  pressure  in¬ 
stantaneously  decays  to  below  the  ambient  pressure.  Therefore, 
the  impulse  is  completely  transmitted  to  the  free-standing  plate 
before  it  recedes  and  the  FSI  has  no  effect  in  reducing  the  blast 
wave  impulse.  Flence,  the  ratio  of  transmitted  impulse  to  the  in¬ 
cident  impulse  lpHj  for  an  infinitesimal  impulse  is  the  same  as 
that  on  a  fixed  wall.  For  a  uniform  incident  blast  wave  with  blast 
intensity  of  50,  the  ratio  of  impulse  on  a  fixed  wall  is  7.2.  The 
reflection  coefficient  is  identical  to  the  ratio  of  impulse  for  a  uni- 


Fig.  9  Comparison  between  the  ratio  of  impulse  /,,//,■  with  and 
without  resistance  for  typical  blast  wave  at  different  plate 
thicknesses 


form  blast  wave.  Since  the  reflection  coefficient  for  a  fixed  wall 
decreases  with  the  decrease  of  blast  intensity,  it  can  be  expected 
that  the  ratio  of  impulse  of  an  exponential  blast  wave  will  be  less 
than  that  of  a  uniform  blast  wave.  It  can  be  seen  from  Fig.  10  that 
the  ratio  of  impulse  for  an  exponential  blast  wave  with  infinitesi¬ 
mal  impulse  is  slightly  below  7.2.  As  the  incident  impulse  in¬ 
creases,  the  free-standing  plate  will  recede  at  a  higher  velocity. 
The  FSI  between  the  blast  wave  and  the  free-standing  plate  results 
in  the  decrease  of  transmitted  impulse  as  well  as  the  ratio  of 
impulse. 

For  incident  peak  overpressures  of  16  MPa  and  32  MPa  and 
decay  time  constants  ranging  from  0.1  ms  to  1.0  ms,  different 
stand-off  distances  are  used  to  obtain  peak  blast  intensities  of  100 
and  200  right  in  front  of  the  free-standing  plate,  where  decay  time 
constants  ranges  from  0.15  ms  to  2.0  ms.  The  ratio  of  impulse  for 
peak  blast  intensities  of  100  and  200  behaves  the  same  as  that  for 
peak  blast  intensity  of  50.  It  can  be  seen  from  Fig.  10  that  in¬ 
creased  peak  blast  intensity  results  in  an  increasing  ratio  of  im¬ 
pulse. 

5  Conclusions 

A  numerical  model  of  the  FSI  between  a  blast  wave  and  a 
free-standing  plate  has  been  developed.  Both  the  blast  wave  re¬ 
flection  in  front  of  the  plate  and  the  shock  wave  induced  at  the 
back  of  the  plate  are  considered  in  the  model.  The  numerical 
model  is  validated  by  comparing  it  with  analytical  and  numerical 
results  of  special  cases.  Two  types  of  blast  wave,  uniform  and 
exponential  blast  wave,  are  investigated  in  this  paper.  For  a  uni¬ 
form  blast  wave,  the  FSI  of  a  heavy  plate  resembles  the  blast 
wave  reflection  of  a  fixed  wall.  It  was  shown  that  the  blast  wave 
reflection  of  a  free-standing  plate  of  infinitesimal  mass  results  in 
an  acoustic  reflection  of  a  blast  wave.  The  impact  on  the  plate  is 
identical  to  the  incident  blast  wave  pressure.  The  FSI  of  a  free¬ 
standing  plate  with  finite  mass  falls  between  these  two  limiting 
cases.  As  the  blast  intensity  increases,  the  reflection  coefficient  of 
the  free-standing  plate  asymptotically  approaches  1. 

For  an  exponential  blast  wave,  numerical  results  show  that  the 
effectiveness  of  FSI  in  blast  wave  mitigation  increases  as  the 
thickness  of  the  plate  decreases.  However,  the  effectiveness 
reaches  a  plateau  after  the  thickness  decreases  beyond  a  critical 
value.  This  is  due  to  the  increased  resistance  at  the  back  of  the 
plate.  As  the  thickness  of  the  plate  decreases,  the  impact  of  the 
blast  wave  drives  the  plate  to  ever  higher  receding  velocity,  which 
induces  a  stronger  shock  wave  at  the  back  of  the  plate.  When  the 
shock  wave  strength  at  the  back  of  the  plate  becomes  comparable 
to  the  incident  blast  intensity,  the  effectiveness  of  FSI  in  reducing 
the  blast  wave  impulse  diminishes.  The  model  shows  that  the 
resistance  at  the  back  of  the  plate  plays  an  important  role  in  the 
FSI,  especially  for  relatively  light  plates.  Neglecting  the  resistance 
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in  the  previous  model  [1]  significantly  overpredicts  the  effective¬ 
ness  of  FSI.  Numerical  investigation  of  the  FSI  between  a  blast 
wave  and  a  free-standing  plate  also  reveals  that  the  ratio  of  im¬ 
pulse  is  also  highly  dependent  on  the  incident  blast  intensity  and 
impulse.  For  a  particular  incident  blast  intensity,  a  small  decrease 
in  the  incident  impulse  can  significantly  increase  the  ratio  of  im¬ 
pulse.  For  a  constant  incident  impulse,  the  increased  blast  inten¬ 
sity  results  in  an  increase  of  ratio  of  impulse. 
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Nomenclature 

Roman 

C  —  reflection  coefficient 
F  =  flux  vector 
I  =  impulse  per  unit  area 
M  =  Mach  number 
R  =  gas  constant 
T  =  temperature 
U  =  solution  vector 
W  =  blast  wave  speed 
a  =  speed  of  sound 
e  =  internal  energy  per  unit  mass 
e,  —  total  energy  per  unit  mass 
h  =  thickness 
m  =  mass 
p  —  pressure 
t  —  time 
u  =  velocity 
x  =  Cartesian  coordinate 
Ap  =  pressure  difference 

Greek 

y  =  ratio  of  specific  heats 
p  =  density 

Subscripts 

0  =  initial  condition,  ambient  condition 
b  =  property  at  the  back  of  plate 
g  =  property  of  induced  shock  wave 
i  =  property  of  incident  blast  wave 


p  =  property  of  plate 
r  =  property  of  reflected  blast  wave 
s  =  peak  value 

u  =  property  of  uniform  blast  wave 
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